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Abstract 

Predicting the statistics of realistic wave-chaotic scattering systems requires, in addition to random 
matrix theory, introduction of system-specific information. This paper investigates experimentally 
one aspect of system-specific behavior, namely the effects of short ray trajectories in wave-chaotic 
systems open to outside scattering channels. In particular, we consider ray trajectories of limited 
length that enter a scattering region through a channel (port) and subsequently exit through a 
channel (port). We show that a suitably averaged value of the impedance can be computed from 
these trajectories and that this can improve the ability to describe the statistical properties of the 
scattering systems. We illustrate and test these points through experiments on a realistic two-port 
microwave scattering billiard. 
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PACS Numbers 

05.45.Mt Quantum chaos; semiclassical methods 
03.65.Nk Scattering theory 
03.65.Sq Semiclassical theories and applications 
05.60.Gg Quantum transport 



I. INTRODUCTION 

Random matrix theory (RMT) has achieved substan- 
tial success at predicting short wavelength statistical 
properties of spectra, eigenfunctions, scattering matrices, 
impedance matrices, and conductance of wave-chaotic 
systems (H-lll- By wave-chaotic systems, we mean that 
the behavior of the wave system in the small wavelength 
limit is described by ray orbit trajectories that are chaotic 
0. In practice, the experimental applicability of RMT 
requires consideration of nonuniversal effects. For exam- 
ple, in the particular case of scattering, the scattering 
properties of an open system depend on the coupling 
between the field within the scattering region and the 
asymptotic incoming and outgoing waves connecting the 
exterior to the scatterer. For this problem, researchers 
have developed methods to incorporate nonuniversal cou- 
pling and port-specific effects into the analysis of short 
wavelength scattering data for systems whose closed clas- 
sical counterparts are ray-chaotic [cl-fT^. 

Previous comparisons of experimental data to RMT 
have often employed ensembles of realizations of the sys- 
tem to compile statistics. To create such ensembles, re- 
searchers have typically varied the geometrical configura- 
tion of the scattering region and/or taken measurements 
at several different wavelengths [l2|, 

MM- These 

variations aim to create a set of systems in which none 
of the nonuniversal system details are reproduced from 
one realization to another, except for the effects of the 



port details. Thus, by suitably accounting for the port 
details, it was hoped that only universal RMT proper- 
ties remained in the ensemble data. However, there can 
be problems in practice. For example, in the case of 
geometrical configuration variation, researchers typically 
move perturbing objects inside a ray-chaotic enclosure 
with fixed shape and size [1, [l^, [13 ; or move one wall 
of that enclosure [15|], to create an ensemble of systems 
with varying details. The problem is that certain walls 
or other scattering objects of the enclosure remain fixed 
throughout the ensemble. Therefore, there may exist rel- 
evant ray trajectories that remain unchanged in many 
or all realizations of the ensemble. We term such ray 
trajectories that leave a port and soon return to it (or 
another port) before ergodically sampling the enclosure 
"short ray trajectories" . 

A similar problem arises for wavelength variation in 
which a band of wavelengths is used. Within a wave- 
length range, if a ray trajectory length is too short, then 
the variation of the phase accumulated by a wave fol- 
lowing that trajectory may not be large enough to be 
considered random. In such a case the effects of specific 
(hence, nonuniversal) short ray trajectories will survive 
the ensemble averaging processes described above. These 
problems will make systematic, nonuniversal contribu- 
tions to the ensemble data, and thus consideration of 
short ray trajectories arises naturally in the semiclassical 
approach to quantum scattering theory d, [l^, [11, IT6|423j . 
Such short-ray-trajectory contributions have been noted 
before in microwave billiards fl2 . [isl . [23 | or for quantum 



transport in chaotic cavities id . 



14j and have either con- 



strained or frustrated previous tests of RMT predictions. 

By a "short ray trajectory" we mean one whose length 
is not much longer than several times the characteristic 
size of the scattering region, and which enters the scatter- 
ing region from a port, bounces (perhaps several times) 
within the scattering region, and then returns to a port. 
A "port" is the region in which there is a connection from 
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the scatterer to the outside world. For illustrative pur- 
poses, in what follows we consider systems with either 
one or two ports; as discussed elsewhere, generalization 
of our results to more ports is straightforward [23 - l27| . 
Note that the short ray trajectories we refer to are dif- 
ferent from periodic orbits [13, HE , which are closed 
classical trajectories; short ray trajectories, as defined 
here, are only important for open systems. Although our 
explicit considerations in this paper are for billiard sys- 
tems (i.e., scattering regions that are homogeneous with 
perfectly reflecting walls), presumably these effects may 
also be present with continuous potentials and are not 
limited to billiards. 

In this paper we go beyond the previous treatments 
to explicitly include additional nonuniversal effects due 
to short ray trajectories. Previous work has examined 
short ray trajectories in cases where the system and the 
ports can be treated in the semiclassical approximation 
jisl [2^ or considered the effect on eigenfunction correla- 
tions due to short ray trajectories associated with nearby 
walls [13, HH. The effects of short ray trajectories on 
wave scattering properties of chaotic systems have been 
explicitly calculated before in the case of quantum graphs 
[33 | and for two dimensional billiards Moreover, mi- 
crowave billiard experimental work has extracted a mea- 
sure of the microwave power that is emitted at a certain 
point in the billiard and returns to the same point af- 
ter following all possible classical trajectories of a given 
length [l3]. The Poisson Kernel approach can also be 
generalized to include short ray trajectories [s^ through 
measurement of a (statistical) optical S-matrix. None of 
this prior work developed a general first-principles deter- 
ministic approach to experimentally analyzing the effect 
of short ray trajectories, as we do here. This paper ex- 
pands on a brief report of our preliminary results fs^. 
Our previous work |34l | demonstrated the effect of short 
trajectories in a one-port wave-chaotic cavity, and we 
generalize the short-trajectory correction to a two-port 
system in this paper. Therefore, this generalized correc- 
tion can be applied to multiple port cases by similar pro- 
cedures. Other new ingredients include detailed exami- 
nation of the effect of individual short ray trajectories, 
and extension to the statistics of other wave scattering 
properties. 

The outline of this paper is as follows. In Sec. |lTl we 
summarize a theoretical approach for removing short ray 
trajectories from single-realization data and ensemble- 
averaged data of any wave properties of a wave scattering 
system. A detailed derivation of the results in Sec. Ullis 
given elsewhere (27j . In Sec. IIIIl we describe experiments 
testing the theoretical approach, and we compare the ex- 
perimental results and the theory in Sec. IIVI These com- 
parisons show that compensation for the effects of short 
ray trajectories improves the agreement between RMT- 
based predictions and measured statistical properties of 
ensemble-averaged data and single-realization data. 

While much of our discussion above has emphasized 
the goal of uncovering RMT statistics, we also wish to 



emphasize that many of our results are also of interest 
independent of that goal. In particular, the key step 
in our approach for uncovering RMT statistics is that 
of obtaining a short-trajectory-based prediction for the 
average impedance over frequency and/or configurations 
(Sees. IIVCI and IIVDI) . and we note that this quantity 
is both experimentally accessible and of interest in its 
own right. As another example, in Sec. IIV Al we consider 
the modification of the free-space radiation impedance 
arising in configurations when there are only a few pos- 
sible short ray trajectories, and this situation applies di- 
rectly to many cases where, although reflecting objects 
are present, there is also substantial coupling to the out- 
side. 



II. REVIEW OF THEORY 

In this paper we consider for concreteness an effectively 
two dimensional microwave cavity (described in Sec. IIII[) 
as the scattering system, and this cavity is connected to 
the outside world via one or two single-mode transmis- 
sion lines terminated by antennas inside the scattering re- 
gion, acting as ports. The results and techniques should 
carry over to other physically different wave systems (e.g., 
quantum waves, acoustic waves, etc.) of higher dimen- 
sion and an arbitrary number of ports [27j . Ray trajec- 
tories within the cavity are chaotic due to the shape of 
the cavity. 

Previously most researchers focused on the scattering 
matrix S [3, H, IH-Hlj , which specifies the linear relation- 
ship between reflected and incident wave amplitudes in 
the channels. We focus on an equivalent quantity, the 
impedance Z — Zo{l + S)/{1 — S), because nonuniversal 
contributions manifest themselves in Z as simple addi- 
tive corrections [42]- Here the diagonal matrix Zq gives 
the characteristic impedance of the scattering channels, 
and is measured at chosen reference planes on the trans- 
mission lines as the ratio of the (complex) transmission 
line voltage V to current / (time variation exp{—iujt) 
is assumed). Impedance is a meaningful concept for all 
scattering wave systems. In linear electromagnetic sys- 
tems, it is defined via the phasor generalization of Ohm's 
law a.s V — ZI, and in the case of N channels connected 
to the scatterer, Z is an x A'^ matrix. A quantum- 
mechanical quantity corresponding to the impedance is 
the reaction matrix, which is often denoted in the liter- 
ature as a: = -iZ [l|, H m, m-g^. In what follows our 
discussion will use language appropriate to the electro- 
magnetic context and scattering from a microwave cavity 
excited by small antennas fed by transmission lines (the 
setting for our experiments). 

Our goal is to describe the universal RMT aspects of 
measurement of this impedance matrix (alternatively of 
the scattering matrix). However, raw measurement of 
the impedance have nonuniversal properties because the 
waves are coupled to the cavity through the specific ge- 
ometry of the junction between the transmission lines 



3 



and the cavity. In prior work, the nonuniversal cou- 
pUng effects were parameterized in terms of corrections 
to the impedance Z. The "perfectly coupled" normalized 
impedanc e i£ o removes nonuniversal features due to the 
ports [T2I. ITsj. 

*Co = i?fl'^'(^-«^K)i?fl'^', (1) 

where Z^, = Rr + iXji is the radiation impedance, em- 
bodying the nonuniversal port properties. Rr is the real 
part of Zr, and Xr is the imaginary part. By "perfectly 
coupled" we mean that waves impinging on the cavity are 
fully transmitted to the cavity with no prompt reflection. 
Note that these quantities are N x N matrices for a sys- 
tem with TV ports. Zr is the impedance measured at 
the previously mentioned reference planes on the trans- 
mission lines when the cavity walls are removed, so that 
the outgoing waves launched at the ports never return 
to these ports. Eq. ([T|) removes the effect of the ports 
in the impedance, so the "perfectly coupled" normalized 
impedance i^o of the cavity is all that remains. We note 
that Zr is experimentally accessible through a determin- 
istic (i.e., non-statistical) measurement described in Sec. 

m 

In Hart et al. [27|, it has been proposed that nonuni- 
versal effects due to the presence of short ray trajectories 
can be removed by appropriate modification of the ra- 
diation impedance Zr Zavg = Zj^-I- [short trajectory 
terms], where Zavg = Zq{1 + — S) stands for the 

"average" impedance. S is the average scattering matrix, 
and the Poisson kernel characterizes the statistical distri- 
bution of the scattering matrix S in terms of S 0, [HI, ■ 
More specifically, the averaged impedance is written as 

Zavg ~ Zr + R](^zR]1^, (2) 

where z is the short-trajectory correction. With 
this modification, i^o in Eq. ([T|) is extended to a 
perfectly-coupled and short-trajectory-corrected normal- 
ized impedance 

- KvT - ^^-s) RavT ' (3) 

where Zavg = Ravg + 'iXavg- Wc define z = p + ix, so 
Ravg and Xavg cau be written as 

Ravg = Rr. + Rr^RRr^j (4) 
Xavg — Xr + R](^ xRr^ ■ (5) 

In the lossless case p and x are the real and imaginary 
parts of z; Ravg and Xavg are the real and imaginary 
parts of Zavg- However, with uniform loss (e.g., due to 
an imaginary part of a homogeneous dielectric constant 
in a microwave cavity), Ravg and Xavg are the analytic 
continuations of the real and imaginary parts of the loss- 
less Zavg as fc — )■ k + ik/{2Q), where Q 3> 1 is the quality 
factor of the closed system, and k denotes the wavenum- 
ber of a plane wave. These analytic continuations are 



no longer purely real (i.e., p, Xj Ravg and Xavg become 
complex). The normalized impedance of the lossy sys- 
tem has a universal distribution which is dependent only 
on the ratio fc/(2QAfc), where Ak is the mean spacing 
between modes [1, [ll, [ll, [13, H SI] . 

For the system with TV ports, the elements of the ma- 
trices z are [13] 



Zn,rn — ^ {-P&(n,m) y -^^(n.m) GXp [- (ifc + K)L6(„,m) 
6(n,m) 

— ikLpQrt{n,m) ~ ^Pb{n,m)^]}- (6) 

With pn,m + iXn,m = Zn,m, the analytic continuations of 
Pn,m and Xn,7n are 

Pn,m = ^2 {~Pb{n,m)\J Di,(^n,m) COs[-(fc — in)Li,^n,m) 
b(n,m) 

— kLpQj.f(^„^jj-^-^ — I3b{n,m)^]}j (7) 



Xn,m = ^ {-Pb{n,m)y Db(^n,m) sin[-(fc - in)L^n,m) 
b{n,ra) 

— kLport(n,m) ~ /^b(n,m) """l }: (8) 

where 6(n, m) is an index over all classical trajectories 
which leave the nth port, bounce Pb(n.m) times, and re- 
turn to the mth port, L^n,m) is the length of the trajec- 
tory b{n,m), K = k/{2Q) is the effective attenuation pa- 
rameter taking account of loss, and Lpart(n,m) is the port- 
dependent constant length between the nth port and the 
mth port. Di,(^n,m) is a geometrical factor of the trajec- 
tory, and it takes into account the spreading of the ray 
tube along its path. This geometrical factor is a function 
of the length of each segment of the trajectory, the angle 
of incidence of each bounce, and the radius of curvature 
of each wall encountered in that trajectory; it has been 
assumed that the port radiates isotropically from a lo- 
cation far from the two-dimensional cavity boundaries. 
Pb{n,m) is the survival probability of the trajectory in the 
ensemble and will be discussed subsequently. In the loss- 
less case (k — 0), note that p and x are both real, but 
they are both complex in the presence of loss. In arriving 
at Eqs. (O and ([5]) it is assumed that the lateral walls 
present perfect-metal boundary conditions, and that foci 
and caustics are absent. These assumptions are well sat- 
isfied for the cavity shape used in our experiments. 

Note that the sums in Eqs. ([7|) and ([S]) involve terms of 
the form sine and cosine of [(fc— iK)Lb(n,m) + - • ■] which for 
large L^n,m) increase exponentially like exp(«;Lb(n,m))- 
Although in the experiments we have observed that the 
parameters Pb{n,m) and Di,(^n,m) decrease exponentially 
when i{,(„^r„) increases, the sums do not necessarily con- 
verge if the loss is too high. Accordingly, we will use 
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a finite cutoff of the sum and regard the cutoff result 
as being asymptotic. In contrast, if instead of extract- 
ing RMT statistics, we regard our goal as approximat- 
ing Zavg ~ Ravg + iXavg, then Eq. ^ shows that we 
only need to consider Zn^m- Thus the sum involved in 
the calculation of z (Eq. ([6])) is now over terms that 
decrease exponentially with increasing path length as 
exp(— kL{,(„ „i)). This sum is much more likely to con- 
verge than the sums in Eqs. ([7]) and ((8)). 

In practice, when considering either of the sums in 
Eqs. (|ni), © or ([5]), we employ a cutoff by replacing the 

sums by which signifies that the sum is now over 

b(n,m) 

all trajectories b{n, m) with lengths up to some maxi- 
mum length Lm, -^^b(n,m) ^ Lm, and this is the rea- 
son of the name "short- trajectory correction". We use 

JLm) ^{Lm) JLm) nC{LM) vi^M) T3(Lm) 1 y(Lm) 

P , X I ^ I I ^avg , J^avg , anu ^i-avg 

to indicate the finite length versions of those quantities. 
Combined with the measured radiation impedance, these 
corrections can be analytically determined by using the 
ray-optics of short ray trajectories between the ports and 
the fixed walls of the cavity. After understanding the 
nonuniversal effects of radiation impedance and short tra- 
jectories, we compare the predictions of RMT with the 
normalized impedance which contains the remaining 
features of longer trajectories and the deviations between 
a single realization and the ensemble average. 

In Sec. mil in addition to configuration averaging 
[taken into account by the quantity Pb(n.m) Eqs. ([B]), 
([7]) or ([5]) (see Sec. IIII|)]. we will also employ frequency 
smoothing. In particular, if /(w) denotes a frequency de- 
pendent quantity, then we take its frequency smoothed 
counterpart to be the convolution of /(w) with a gaus- 
sian, 

7H = y - ^')duj', (9) 

where 

Applying the operation ([9]) to our short-trajectory correc- 
tion formulas, Eqs. dU, ^ or ([8]), with k ^ uj/c, we see 
that the summations acquire an additional multiplicative 
factor, 

exp[-^(ib(«.m) + ^port(n,m))^(Aw/c)^]. (11) 

Thus, as ib(„,„i) increases (i.e., longer trajectories are 
included), the factor ([TT|) eventually becomes small, thus 
providing a natural cutoff to the summations in Eqs. (O, 
dH) and JS]). 

III. EXPERIMENT 

We have carried out experimental tests of the short- 
trajectory effects using a quasi-two-dimensional mi- 



crowave cavity with two different port configurations, a 
one-port system and a two-port system. For the one-port 
case Z and S are scalars, while they are 2x2 matrices in 
the two-port case. Microwaves are injected through each 
antenna attached to a coaxial transmission line of char- 
acteristic impedance Zq, and the antennas are inserted 
into the cavity through small holes (diameters about 1 
mm) in the lid, similar to previous setups [6, 25, 26, 47]. 
The waves introduced are quasi-two-dimensional (cavity 
height 0.8 cm) for frequencies (from 6 to 18 GHz) be- 
low the cutoff frequency for higher order modes (~ 19 
GHz), including about 1070 eigenmodes of the closed 
cavity. The quasi-two-dimensional eigenmodes of the 
closed system are described by the Helmholtz equation 
for the single non-zero component of electric field (Ez), 
and these solutions can be mapped onto solutions of the 
Schrodinger equation for an infinite square well potential 
of the same shape 0, HI] . 

The shape of the cavity walls is a symmetry-reduced 
"bow-tie billiard" made up of two straight walls and two 
circular dispersing walls [48], as shown in the insets of 
Fig. [TJ The lengths of wall A and wall B are 21.6 cm 
and 43.2 cm, respectively, while the wavelength of the 
microwave signals ranges from 1.7 cm to 5.0 cm, putting 
this billiard system just into the semiclassical regime. De- 
spite the fact that we are not deep inside the semiclassical 
regime, we obtain results in good agreement with theory 
(Sec. IIV[) . In the one-port experiments the location of 
the port is 18.0 cm from wall A and 15.5 cm from wall 
B. In the two-port experiments we add the other port at 
the location 35.6 cm from wall A and 15.6 cm from wall 
B. The cavity shape yields chaotic ray trajectories and 
has been previously used to examine eigenvalue [4^ and 
eigenfunction [s^l statistics of the closed system in the 
crossover from GOE to GUE statistics as time-reversal 
invariance is broken. In addition, this cavity was also 
used to study scattering (S-), 0,1431 impedance (^-) [1| 
and conductance (G-) [20] statistics. 

We use an Agilent PNA network analyzer to mea- 
sure the frequency dependence of the complex scatter- 
ing matrix S, and compute the corresponding impedance 
Z = Zo{l + S)/{1 - S). The walls are fixed relative to 
the ports in all experiments. We experimentally deter- 
mine the radiation impedance Zf; by placing microwave 
absorbers along all the side-walls of the cavity and mea- 
suring the resulting impedance. The absorbers eliminate 
reflections from the walls, so this method removes the 
effect of ray trajectories, leaving only the effects of the 
port details. To verify the theory, in preliminary experi- 
ments microwave absorbers are placed only along specific 
walls, and we measure and examine specific contributions 
of individual walls, or groups of walls, to the impedance 
(Sec. IIV Al and lIVBp . In another set of experiments, we 
remove all the absorbers and make measurement on the 
empty cavity (Sec. lIVCp . Finally, we also make mea- 
surement with two perturbers added to the interior of 
the cavity (Sec. lIVPp . Each perturber is a cylindrical 
piece of metal with a height similar to that of the cavity. 
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and their function is to block ray trajectories. The two 
perturbers are systematicahy moved to create 100 dif- 
ferent reaUzations used to form an ensemble of systems 
in which the empty-cavity short trajectories are partially 
destroyed. In the one-port experiments the cross sections 
of the two cylindrical perturbers are irregularly-starlike 
shapes with the maximum diameters 7.9 and 9.5 cm. In 
the two-port experiments the cross sections are identical 
circular shapes with diameter 5.1 cm. 



IV. RESULTS 
A. Individual short ray trajectories 

The first experiment tests whether the theory of short- 
trajectory corrections can predict the effect of individ- 
ual ray trajectories, as well as the aggregate effect of a 
small number of ray trajectories, on the impedance. To 
isolate individual trajectories, microwave absorbers are 
employed to cover some of the cavity walls, and the ex- 
periment systematically includes short ray trajectories 
involving bounces from exposed walls of the 1 /4-bow-tie 
cavity (see insets of Fig. [1] for examples), with no per- 
turbers present. 

Fig. [1] shows comparisons between the measured 
impedance (Zw) and the theoretical form (Z^**-*). Here 
the impedance Z\y is measured from the microwave cav- 
ity with specific walls (W) exposed, where W = B, CD, 
or BC stands for one or more of the walls A, B, C, and D 
shown in the insets of Fig. [TJ We have examined different 
combinations of exposed walls and choose some represen- 
tative cases to show here. Figs. [1] (a) and (b) are for the 
cases of one-port experiments with (a) one wall exposed 
(wall B) and (b) two walls exposed (walls C and D), and 
Figs. [1] (c) and (d) are for cases of two-port experiments 
with (c) one wall exposed (wall B) and (d) two walls 
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exposed (walls B and C). Notice that Z, 
in this case because there is only one single realization. 
Thus, there is no ensemble averaging and Pb(n,m) = 1 for 
all trajectories in Eqs. (O and ([8]). 

Here we focus on the effects of ray trajectories, so we 
remove the effect of the port mismatch from the mea- 
sured impedance. We term this quantity the impedance 
correction 



'1/2 



(12) 

and in Fig. [T] the measured data (red and solid) are 
constructed as the resistance correction Re[zcor] and the 



reactance correction Imfzc 



In each case, the radiation 



impedance of the antenna {Zjf) is determined by a sepa- 
rate measurement in which all four walls are covered by 
the microwave absorbers [6]. According to Eq. ([2]), the 
corresponding theoretical term is 



^cor — ^B. 



(Lm) 
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FIG. 1: (Color online) Plot of corrections to the impedance 
due to (a) one-port with one wall (wall B), (b) one-port with 
two walls (walls C and D), (c) two-port with one wall (wall B), 
and (d) two-port with two walls (walls B and C) short trajec- 
tories, versus frequency from 6 to 10 GHz inside the 1/4-bow- 
tie cavity. Shown are the theoretical predictions (Re[2:co r'*^'] 
and Im[2:cor '^'l) in blue (dashed) curves and experimental data 
(Re[2:cor] and Im[2:cor]) in red (solid) curves. 
The insets show the 1/4-bow-tie billiard with thicker lines 
representing walls covered by microwave absorbers, the ports 
are shown as circular dots, and a few short trajectories are 
represented with lines as illustrations. 
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Therefore, the theoretical curves (blue dashed) of the re- 
sistance correction and the reactance correction in Fig. [T] 

are Re[zio^^] and Ini[2;cor"^], respectively, and they can 
be calculated by using Eq. with the known geometry 
of the cavity and port locations, including all short orbits 
up to the maximum length Lm = 200 cm. 



The measured data generally follow the theoretical pre- 
dictions quite well, thus verifying that the theory offers 
a quantitative prediction of short-trajectory features of 
the impedance Zw In Figs. [1] (a) and (b), the one- 
port cases, we examine the resistance (reactance) correc- 
tion for the waves entering and returning from the cavity 
through the single port, and the short-trajectory correc- 
tions respectively include one trajectory ton W = B and 
a sum over nine trajectories for W = CD. In Figs. [T] (c) 
and (d), the two-port cases, we examine the resistance 
(reactance) correction between the two ports. It corre- 
sponds to the elements Zcor, 1,2 and in the 2x2 

matrices (zcor and Zcor'^), and the short-trajectory cor- 
rections respectively include sums over two trajectories 
(W — B) and four trajectories {W = BC). For ihustra- 
tion, the insets of Fig. [T]show some representative short 
ray trajectories. Note the direct trajectory between the 
two ports without bouncing on walls is treated as a ray 
trajectory in Eqs. ([7]) and 



For the propagation attenuation, a frequency- 
dependent attenuation parameter k(/) is calculated uti- 
lizing the previously measured frequency-dependent loss 
parameter a [5l| for this cavity. Fig. [2] shows the 
loss parameter a{f) and the attenuation parameter k(/) 
for the empty cavity case. Here a is defined as the 
ratio of the closed-cavity mode resonance 3-dB band- 
width to the mean spacing between cavity modes, a — 
P/{Ak^Q) ~ PA/{4nQ), where A is the area of the 
cavity {A = O.llSm^), k is the wavenumber, and Q is 
the quality factor of the cavity. The mean spacing be- 
tween modes varies from 21 MHz at 6 GHz to 6.9 MHz at 
18 GHz. We obtain the loss parameter a from the mea- 
sured impedance data in the empty cavity case according 
to the procedures presented in previous work IH, [5l| . 
In this procedure the frequency-dependent loss parame- 
ter a{f) is determined by selecting a proper frequency 
range (e.g., 1.8 GHz), computing the PDF of the "per- 
fectly coupled" normalized impedance i^o (Eq. ((Ij), and 
then comparing to different PDFs that were generated 
from numerical RMT, using a as a fitting parameter. 
Once a{f) is known, the frequency-dependent attenua- 
tion parameter k(/) can be calculated because the domi- 
nant attenuation comes from losses in the top and bottom 
plates of the cavity, and it is well modeled by assuming 
that the waves suffer a spatially uniform propagation loss 
AC = fc/(2Q) = 2TTa/{kA). 




0.08-1 ■ . ■ , ■ , ■ , ■ . ■ . 
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FIG. 2: (Color online) Plot of the loss parameter a and the 
attenuation parameter n versus frequency inside the empty 
1/4-bow-tie cavity. 



B. Sources of errors 

We propose that there are two major sources of the 
deviations between the theory and experiment shown in 
Fig. [T] The first is that the microwave absorbers do 
not fully suppress the trajectories, and the second arises 
from the ends of microwave absorbers that scatter energy 
back to the ports. To verify this, the k — lo/c dependent 
impedance corrections data in Fig. [T] are Fourier trans- 
formed to the length domain and shown in Fig. [3] as 

e{L)^eas = c\FT{z,or}\ and e{L)theory = c|FT{z^^/^)}|, 
where c is the speed of light, and FT{...} is the Fourier 
transformation (k ^ L). Note that the frequency range 
of the Fourier transformation is from 6 GHz to 18 GHz, 
and therefore, the resolution in length is 2.5 cm. 

In the length domain, the major peaks of the mea- 
sured data (red) match the peaks of the theoretical pre- 
diction (blue and dashed), and this verifies that the short- 
trajectory correction can describe the major features of 
the measured impedance in the scattering system. For 
example, the matched peak in the theory curve and the 
data curve in Fig. [3Ka) corresponds to the short trajec- 
tory from the port to wall B and returning, shown as the 
red (vertical) line in the inset of[3](a). However, there are 
several minor peaks in the measured data not present 
in the theoretical curves. After further examination of 
the geometry, the positions of these deviations in Fig. [3] 
match the lengths of trajectories which are related to the 
partially-blocked corners of the cavity, or bounce off mi- 
crowave absorbers with a large incident angle. When the 
microwave absorbers end at the corners, they produce 
gaps and edges, and these defects create weak diffractive 
short trajectories. For example, the green lines El and 
E2 in the insets in Figs. [Sja) and (b) and E2 in Figs. [3jc) 
and (d) represent the diffractive short trajectories leav- 
ing a port, bouncing off the partially covered corners, 
and returning to a port. Their path lengths match the 
deviations between the measured data and the theory as 
labeled in the figure. Furthermore, the blue line E3 in 
the inset in Fig. ^ih) represents the short trajectories 
produced by the corner and bounced from one wall. 
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FIG. 3: (Color online) Plot of the absolute values of the 
impedance correction in the length domain e{L) in the cases 
of (a) one-port with one wall (wall B), (b) one-port with two 
walls (walls C and D), (c) two-port with one wall (wall B), 
and (d) two-port with two walls (walls B and C) short tra- 
jectories, versus length from to 120 cm inside the 1/4-bow- 
tie cavity. Shown are the measured data in red curves with 
circular symbols and the theory in blue dashed curves with 
triangular symbols. The insets show the 1/4-bow-tie billiard 
with labeled trajectories that show the sources of error (El, 
E2, and E3). 



The other source of error is imperfection of the mi- 
crowave absorbers that reflect ^ — 20 dB of the incident 
signal for normal incidence, and more for oblique inci- 
dence. Therefore the short trajectories shown as purple 
hues El in Figs. [3jc) and (d) with a large incident angle 
bring about deviations between the theory and exper- 
imental data. These sources of error due to the ends 
of absorbers or large incident angles on absorbers were 
not included in the short-trajectory correction. However, 
these errors will not concern us further because all mi- 
crowave absorbers are removed from the cavity in the 
experiments discussed below. 

Besides the errors discussed above, another source of 
error is the difficulty in reproducing the antenna ge- 
ometry with each measurement as the cavity is opened 
and re-sealed between the measurement of the radiation 
impedance and exposed wall cases. Another concern is 
multiply-reflected trajectories that bounce off of the an- 
tennas. However, because we describe trajectories in 
terms of the impedance instead of the scattering param- 
eter, the multiply-reflected trajectories are incorporated 
in a single impedance term (see [13] )■ This is an impor- 
tant advantage of using impedance because it can take 
account of the multiple-reflected trajectories in a simple 
compact form. 



C. Short-trajectory correction in the empty cavity 

In the next experiment, an empty 1/4-bow-tie cav- 
ity with no microwave absorbers or perturbers is used. 
Therefore, all possible trajectories between ports are 
present in this single realization. Fig. 2] shows (a) the 
real and (b) the imaginary parts of the first diagonal 
component of the impedance in a two-port cavity, cor- 
responding to Rc[Zii] and Im[Zii]. Figs. S] (c) and 
(d) are for the off-diagonal component of the impedance, 
Re[Zi2] and Im[Zi2]. The radiation impedance (black) 
traces through the center of the fiuctuating impedance 
data of the empty cavity and represents the nonuniversal 
aspects of the coupling antennas @, [H, ■ Note for the 
Zi2 case, all signals from port 1 to port 2 are treated as 
trajectories, so the off-diagonal radiation impedance is 
zero. Also shown in Fig. 21 the theoretical impedance 
Z(Lm) (Eqs. dll) and dH)) include a finite number of 
short trajectories {Lm = 200 cm, given a total of 584 

trajectories for z[^'^''^ and 1088 trajectories for z'^'^''^). 
Note that the attenuation parameter k(/) is determined 
through the same procedures as in Sec. IIV Al 

In Fig. m the theoretical impedance Z^-^*^) tracks 
the main features of the single-realization measured 
impedance although there are many sharp deviations be- 
tween the two sets of curves. These fluctuations are 
expected because of the inflnite number of trajectories 
(-^fc(n,m) > 200 cm) not included in the theory. 

It is more appropriate to compare the theory to a fre- 
quency ensemble of single-realization data. A frequency 
ensemble is created by considering frequency smoothed 
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FIG. 4: (Color online) Plot of the impedance due to the tra- 
jectories of four side-walls, versus frequency from 6 to 8 GHz 
inside the empty 1/4-bow-tie cavity. Shown are (a) the real 
part of Zii, (b) the imaginary part of Zn, (c) the real part 
of Zi2 and (d) the imaginary part of Zi2. The three curves 
are the measured impedance Z (red thinner), the theoreti- 
cal impedance Z^'^'^'^ (blue dashed), as well as the radiation 
impedance Zr (black). 
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FIG. 5: (Color online) Plot of the smoothed impedance versus 
frequency from 6 to 8 GHz. Shown are, (a) for Zn and (b) 
for Zi2 , the real (three upper curves) and the imaginary part 
(three lower curves) of the smoothed impedance for the theory 
(^^{Lm) ^j^jj Lm = 200 cm, blue dashed) and the experiment 
(red solid), as well as the measured (un-smoothed) radiation 
impedance of the port {Zr, black thick). 



experimental data and comparing it vifith the smoothed 
theoretical prediction for system-specific contributions to 
the impedance. The frequency smoothing (Eqs. ^ and 
([To)) ) suppresses the impedance fluctuations due to long 
trajectories and reveals the features associated with short 
trajectories. Fig. [S]shows the radiation impedance (black 
thick), the smoothed measured impedance Z (red solid) 
and the smoothed theoretical impedance Z*^-^"' (blue 
dashed) . The smoothing is made by a Gaussian smooth- 
ing function with the standard deviation Aa;/(27r) = 240 
MHz (Eq. PU]) ). Gaussian frequency smoothing inserts 
an effective low-pass Gaussian filter on the trajectory 
length, and thus, the components of impedances {Z and 
Z(Lm)^ in the length domain are limited to those with the 
path length L < 27rc/Aa; = 125 cm. Fig. [5] shows that 
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the smoothed theory matches the similarly smoothed ex- 
perimental data very well, therefore, the short trajectory 
theory correctly captures the effects of ray trajectories 
out to this length region (< 125 cm). 

When computing the sum of short-trajectory correc- 
tion terms in a low loss case like the empty cavity, a 
problem appears associated with the finite number of 
terms (-/jb(„,m) ^ Lm) in the sum over trajectories. In 
the theory, to perfectly reproduce the measured data in 
the empty cavity requires an infinite number of correc- 
tion terms. Therefore, in some frequency regions where 
the experimental impedance changes rapidly, the finite 
sum for the theoretical resistance Re[Z(^"'^^] will show 
values less than zero, which are not physical for a passive 
system. For example, see the blue dashed curve in Fig. 
mja) between 7.3 and 7.4 GHz. This problem is similar to 
Gibbs phenomenon in which the sum of a finite number 
of terms of the Fourier series has large overshoots near a 
jump discontinuity. 

Besides verifying that the short-trajectory correction 
agrees with system-specific features of the measured data, 
we would next like to demonstrate that including short- 
trajectory corrections improves the ability to reveal un- 
derlying universal statistical properties, even in a single 
realization of the system. We compute the statistical 
properties of the real and imaginary part of the normal- 
ized impedance difference 



R 



-1/2 



-1/2 



(14) 



in 500 MHz frequency windows from 6 to 18 GHz for a 
single realization of the bow-tie cavity. Random matrix 
theory predicts that the distribution of Re[z^^*^''] and 

Im[z^^"^] should have zero means, and identical standard 
deviations d, H, [H, HI, IH, H^l . For the one-port case, we 
use the normalized impedance difference value directly, 
and for the two-port case, we consider the eigenvalues of 
the 2x2 normalized impedance difference matrix. 

Fig. [S] shows that the RMS errors of statistical param- 
eters between the measured data and the theory decrease 
upon including more short trajectories in the correction 
(i.e., increasing Lm)- We compute the root mean square 
value of errors (C defined later) for a series of frequency 
windows covering the range from 6 to 18 GHz, and the 
results are shown versus different short-trajectory cor- 
rections with varied maximum lengths Lm- Notice the 
cases of Lm = denote the impedance corrected by only 
the radiation impedance without any short-trajectory 
correction. For the normalized impedance difference 
z^^"'^"* in each frequency window, the RMS error is 
defined as the root mean square value of — 0| or 
l/ix — 0| for the difference of the measured mean from 
zero, and C,a is defined as the root mean square value 
of \aR — axMipR + Ox) for the difference of standard 
deviations between the real part PDF and the imaginary 
part PDF. Here /i^ and are respectively the means of 
Re[zf ^^^] and \^\zf"'\ in each window, and ub. and ax 
are the standard deviations of Re[z^^"''] and Im[z^^"-'] in 
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FIG. 6: (Color online) Plot of the RMS errors for (a) the one- 
port experiment and (b) the two-port experiment. Shown are 
C^fi for the mean of the real part (red circles) and the imaginary 
part (blue triangles) of the normalized impedance difference 
and ("o- for the difference of standard deviations (green rectan- 
gles), versus short-trajectory corrections with the maximum 
length from to 200 cm. Data are taken for a single realiza- 
tion of the bow-tie cavity. 



each window. In all three statistical parameters for both 
one-port and two-port cases, the RMS errors C decrease 
when we correct the data with more short trajectories. 
This verifies that using short-trajectory corrections in a 
single realization of the wave-chaotic system can more ef- 
fectively reveal the universal statistical properties in the 
data. 



D. Configuration averaging approach 

Many efforts to determine universal RMT statistics in 
experimental systems are based on a configuration av- 
eraging approach that creates an ensemble average from 
realizations with varied configurations. In principle, one 
can recover the nonuniversal properties of the system 
@ El via the configuration averaging approach, which 
is motivated by the "Poisson kernel" theory of Mello, 
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Pereyra, and Seligman fvll. Specifically, ensemble aver- 
ages of the measured cavity data are used to remove the 
system-specific features in each single realization. Note 
that in the past, the configuration averaging approach 
was explicitly assumed to only remove the effects of the 
nonuniversal coupling; however, it was recently general- 
ized to include the nonuniversal contributions of short 
trajectories 



Here the experimental results verify that the short- 
trajectory correction (Eqs. ([7]) and ([8])) can describe 
nonuniversal characteristics of wave-chaotic systems in 
the configuration ensemble. Two cylindrical pieces of 
metal are added as perturbers in the wave-chaotic sys- 
tem that is shown in the inset of Fig. [71 where the 
dots represent the ports and the two circles represent 
the perturbers. The locations of the two perturbers in- 
side the cavity are systematically changed and accurately 
recorded to produce a set of 100 realizations for the en- 
semble [i^ 0, 113] • The scattering matrix S is mea- 
sured from 6 to 18 GHz, covering roughly 1070 modes 
of the closed cavity. After the ensemble average, longer 
ray trajectories have higher probability of being blocked 
by the two perturbers in the 100 realizations; therefore, 
the main nonuniversal contributions are due to shorter 
ray trajectories. We compare the measured ensemble 
averaged impedance (Z) and the theoretical impedance 
that is calculated from Eqs. ^ and ^ with the 



7(Lm) 

^avg 



maximum trajectory length Lm = 200 cm. 

In the configuration ensemble, contrary to the previous 
cases without perturbers, we need to introduce a survival 
probability Pb{n,m) for each ray trajectory term in Eqs. 
([7|) and ([8|) . Notice that the two perturbers can block ray 
trajectories and influence their presence in the ensemble 
realizations. Thus, we multiply each term in the sum 
by a weight Pb{n,m) equal to the fraction of perturbation 
configurations in which the trajectory is not intercepted 
by the perturbers. The values of Pb{n,m) are between 
and 1, and a longer ray trajectory generally has a higher 
chance of being blocked by perturbers, and thus it has 
smaller Pb(n.,m)- 

Note that by recording the positions of perturbers in 
all realizations, we are able to do impedance normaliza- 
tion with short-trajectory correction individually for each 
realization, similar to the procedures in the empty cav- 
ity case. Here we introduce Pb{n,m) as a more general 
description for the case in which only the probabilities 
of survival of particular short trajectories are known. In 
addition, we ignore the effect of newly created trajecto- 
ries by the perturbers in each specific realization because 
they are averaged out in the ensemble. Note that the 
attenuation parameter k in the short-trajectory correc- 
tion terms (Eqs. ([7]) and ([S])) is recalculated using the 
measured PDFs of impedance in the ensemble case, using 
procedures similar to those for the case of the empty cav- 
ity. Due to the presence of two perturbers in the cavity, 
the attenuation parameter and loss parameter are slightly 
larger ('^ 0.1 for a) than in the empty cavity case. 

The result of comparisons between the 2x2 matrices 
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FIG. 7: (Color online) Plot of the average impedance ver- 
sus frequency from 6 to 8 GHz. Shown are, (a) for Zn and 
(b) for Zi2, the real (three upper curves) and the imaginary 
part (three lower curves) of the average impedance for the 
theory {Zivg''' with Lm ~ 200 cm, blue dashed) and the con- 
figuration average experiment {{Z) red solid), as well as the 
measured radiation impedance of the ports {Zr, black thick). 
Plot (c) shows the real part of the measured impedance Z12 
(black thin) in a single realization comparing with Z^J^g'12 and 
{Zi2). Inset: The wave-chaotic two-dimensional cavity with 
perturbers and two ports. 



11 



{Z) and Zavg of the two-port experiment is shown in 
Fig. [T] We have pubhshed the resuh of the one-port 
experiment in another paper [s^l- Here Fig. [7] (a) shows 
the comparison of the first diagonal component of the 
impedance, and Fig. [7] (b) shows the comparison of the 
off-diagonal component. The measured data (red solid) 
foUow the trend of the radiation impedance (black thick) , 
and the theory (blue dashed) reproduces most of the fluc- 
tuations in the data by including only a modest num- 
ber of short-trajectory correction terms. Fig. illus- 
trates the comparison between the measured impedance 
Z in a single realization and (Z) The strong fluctua- 
tions in the measured impedance Z (black thin) curve 
have diminished due to the configuration averaging, and 
the remaining fluctuations of the configuration averaged 
impedance (Z) away from the radiation impedance are 
closely tracked by the theoretical (blue dashed) curve. 
Note that no wavelength averaging is used here. The 
good agreement between the measured data and the the- 
oretical prediction verifies that the new theory, Eqs. ([7]) 
and ([5]), predicts the nonuniversal features embodied in 
the ensemble averaged impedance well. 

The deviations between the (Z) curves and the ziyg'^ 
curves in Fig. [7] may come from several effects. The first 
is the remaining fluctuations in the (Z) due to the finite 
number of realizations. We estimate this to be on the 
order of a^^z) — (Jz/V^OO ~ Ifl, where (7(2) is the stan- 
dard deviation of ensemble averaged impedance (Z) of 
100 realizations, and az is the standard deviation of the 
measured impedance Z of a single realization. This ac- 
counts for the remaining sharp features in (Z). Another 
source of errors is that we ignore the effect of newly cre- 
ated trajectories by the perturbers. Even though these 
new trajectory terms are divided by 100 (the number 
of realizations), the remaining effects create small devia- 
tions. 

We have demonstrated that the ensemble average tech- 
nique reveals the nonuniversal properties of the system, 
and short-trajectory corrections can help to better de- 
scribe the nonuniversal part than using the radiation 
impedance alone. Now we demonstrate the benefit of 
using the short-trajectory correction in another respect. 
That is, when removing the nonuniversal part to reveal 
the universal properties, the revealed universal proper- 
ties have better agreement with the prediction of RMT 
if we take account of the short-trajectory correction. We 
compare the probability distribution of the eigenvalues 
of the normalized impedance i^o [normalized by the ra- 
diation impedance (see Eq. (UJ)], [normalized by 
the short-trajectory-corrected impedance (see Eq. ^)], 
and the PDFs that were generated from numerical RMT 

Fig. |5] shows the PDFs of the normalized impedance 
i^Q, and the corresponding data from numeri- 

cal RMT. For i^o in Eq. ([T]), we take the impedance 
Z from the 100 realizations and in a frequency range of 
200 MHz. Similarly, for in Eq. ©, we use the 

same measured data but consider short-trajectory cor- 
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FIG. 8: (Color online) Plot of probability distributions of (a) 
the real part and (b) the imaginary part of the normalized 
impedance eigenvalues in the frequency range 6.8 GHz ~ 7.0 
GHz, and (c) the real part and (d) the imaginary part of 
the normalized impedance in the frequency range 11.0 GHz 
~ 11.2 GHz. The blue (squares) curve is that is nor- 

malized with the short-trajectory-corrected impedance, the 
green (circles) curve is j^o that is normalized with the radi- 
ation impedance only, and the black curve shows the PDF 
generated from numerical RMT. 



rection with Lm = 200 cm in addition to the radiation 
impedance. Here we show two examples for frequency 
ranges 6.8 GHz - 7.0 GHz in Figs. |H| (a) and (b); and 
11.0 GHz ~ 11.2 GHz in Figs. |5] (c) and (d). It should 
be noted that in the past such narrow frequency ranges 
were not used because short trajectories created strong 
deviations from RMT predictions d, [H, li^l- Because 
the normalized impedance is complex. Figs. [8] (a) and 
(c) show the distribution of the real part of the normal- 
ized impedance; (b) and (d) are the imaginary part. For 
the numerical RMT data, the loss parameters [a — 0.3 
for 6.8 GHz - 7.0 GHz and a = 0.4 for 11.0 GHz - 11.2 
GHz) were determined by the best matched distribution 
with a much wider frequency range (2 GHz). As seen 
in Fig. [5] the distribution of the normalized impedance 
has much better agreement with the prediction of RMT 
when we consider the short-trajectory correction up to 
the trajectory length of 200 cm. 

The deviations of the distribution between the nor- 
malized impedance and the predictions of RMT in Fig. 
[8] are due to short trajectories that remain in the ensem- 
ble average of the 100 realizations. We have seen that 
the PDFs of and approach each other when 

we use the data in a wider frequency range. This is be- 
cause in a wide enough frequency window the fluctuations 
in the impedance due to a trajectory can be compen- 
sated (i.e., the required window is 1.8 (GHz) ~ c/17(cm) 
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for the shortest trajectory with i) = 15 cm and 
Lport{i,i) = 2 cm). Therefore, we take a 2 GHz fre- 
quency window to obtain a universal distribution that is 
independent of which normahzation methods we use (Eq. 
([1]) or Eq. ([3]) ) , and we also determine the loss parameter 
from this universal distribution [2^ [sij . 

Furthermore, the deviations shown in Fig. [5] match 
the difference between (Z) and Zn shown in Fig. [T] For 
example, in the frequency range 6.8 GHz ~ 7.0 GHz the 
ensemble averaged impedance (Zn) is smaller than the 
radiation impedance Zr^h in the real part and larger 
in the imaginary part, and in Figs. [5] (a) and (b) we 
can see the same bias of the distribution of the normal- 
ized impedance i^o that is normalized with the radiation 
impedance only. Therefore, with short-trajectory correc- 
tions, we can better explain the deviations between the 
measured ensemble average and the universal properties 
predicted by RMT. 



E. Uncovering RMT statistics of the scattering 
matrix 



We now test the benefits of the short-trajectory cor- 
rection in uncovering universal RMT statistics of the 
scattering matrix. We find that including short tra- 
jectories in the impedance normalization improves our 
determination of the RMT statistical properties of the 
scattering matrix S, given the same amount of ensemble 
and frequency averaging. In other words, we show that 
the need to resort to wavelength averaging over large 
numbers of modes is significantly reduced after includ- 
ing short-trajectory corrections in the impedance nor- 
malization. This section shows the result of the two-port 
experiment; the result of the one-port experiment has 
been published in our previous work [3J|. Here we ex- 
amine the phase tps of the eigenvalues of the normalized 
scattering matrix s(^") because the statistics of ips do 
not change with the frequency-dependent loss parameter 
c^(/) [111- The normalized scattering matrix s(^Af) is de- 
fined as - l)/(i^(^*^) + 1) = |s(^")|e''^% where 
the normalized impedance is given by Eq. ([3]) with 
corrections by Eqs. ([7]) and dS]). Note for the two-port 
experiment, we compute the eigenvalues of the matrix 



and ifs is the phase of the eigenvalues. RMT pre- 



diets that ips should have a uniform distribution from 
to 2n independent of loss 0, Q i verified in previous 
experiments (25ll47|. 

Fig. HI illustrates the benefits of using short-trajectory- 
corrected data to examine the statistical properties of the 
scattering matrix. An ensemble of data is created by av- 
eraging 100 realizations of the cavity with the two per- 
turbers present, between 6 and 18 GHz, encompassing 
about 1070 closed-cavity modes. Theoretical short tra- 
jectory corrections are weighted by the survival probabil- 
ity iPb{n,m)) which are calculated from the known loca- 
tions of the perturbers. Fig. IHl^a) shows a representative 
probability distribution {P{ips)) of the phase of eigenval- 
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FIG. 9: (Color online) Plot of (a) an example probability 
distribution of the phase ifa of eigenvalues of the normal- 
ized scattering matrix s'^*^' from ~ 2-k, taken over the 
11~11.5 GHz frequency window inside the 1/4-bow-tie cavity 
of 100 realizations with perturbers. The blue (squares) curve 
is from data normalized with Lm ~ 200 cm, the green (circles) 
curve is from data normalized with the radiation impedance 
only, and the black line shows the perfectly uniform distri- 
bution for comparison. Plot of (b) the average RMS error 
of distributions of (fis, where the normalized scattering ma- 
trix s'^m) jg calculated from impedance normalized with the 
radiation impedance only (green circles), with the measured 
ensemble averaged impedance (pink stars), or with theoret- 
ical impedance with Lm = 50 cm (red triangles), 100 cm 
(orange invert triangles), up to 200 cm (blue squares), versus 
frequency window size from 10 MHz to 1.0 GHz. 



ues of the normalized scattering matrix. The data for 
the normalized scattering matrix are formed by taking 
the measured impedance from each of the 100 realiza- 
tions and normalizing with the radiation impedances Zji 
or the theoretical impedance with short-trajectory cor- 
rection and these PDFs are computed from the 
data in a 500 MHz frequency range (11.0 ~ 11.5 GHz). 
In this frequency range, it is clear that the distribution 
normalized with short trajectories is significantly closer 
to a uniform distribution than the one normalized with- 
out short trajectories. 
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To quantify this improvement, it is convenient to look 
at the average root-mean-square (RMS) error with re- 
spect to the uniform distribution for the two types of 
normahzation (Fig. (Hllb)). Here the RMS error is de- 
fined as 



RMS error — 




where rii is the number of elements in the i bin in the 
10 bins histogram of P{(fs), and (n^) is the mean of n^. 
Therefore, when a distribution is more uniform, its RMS 
error is smaller. The average RMS error of each fre- 
quency window is averaged over the spectral range from 
6 to 18 GHz, and Fig. [SJb) shows the average RMS 
error versus the frequency window sizes from 10 MHz 
to 1.0 GHz. This figure shows that the distributions 
of data are systematically more uniform as more short 
trajectories are taken into account in the normalization 
for a given window size. The average RMS error of the 
data normalized without short trajectories (the radia- 
tion impedance only (green circles)) has the largest error. 
When we include more short trajectories from the maxi- 
mum lengths Lm — 50 cm (red triangles) up to Lm — 200 
cm (blue squares), the average RMS errors decrease. The 
improvement is obvious after including short trajectories 
with Lm = 50 cm, and it saturates for trajectories with 
Lm = 100 cm. For comparison, we also add a curve (pink 
stars) for the data normalized by the measured ensemble 
averaged impedance (Z) (the red curves in Fig. [7]), and 
it is the most uniform case. It is observed that the im- 



provement of statistical properties with short-trajectory 
corrections is more significant when the frequency win- 
dow size is smaller. 



V. CONCLUSIONS 

A theory for the nonuniversal effects of coupling and 
short ray trajectories on scattering wave-chaotic systems 
has been tested experimentally on a two-dimensional, 
wave-chaotic, electromagnetic cavity with one or two 
transmission line channels connecting to the outside. In 
particular, the theoretical predictions match the mea- 
sured data in the cases of a frequency ensemble in a single 
realization and the configuration ensemble in the ensem- 
ble of 100 realizations. By removing nonuniversal effects 
from measured data, we can reveal underlying universally 
fluctuating quantities in the scattering and impedance 
matrices. These results should be useful in many fields 
where similar wave phenomena are of interest, such as 
nuclear scattering, atomic physics, quantum transport in 
condensed matter systems, electromagnetics, acoustics, 
geophysics, etc. 
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